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Abstract 

Quarkonium measurements in heavy-ion collisions play an essential role in understanding the hot, dense medium 
created in such collisions. As a reference, their production mechanism in p -l- p collisions needs to be thoroughly 
understood. In this paper, we report the measurement of inclusive cross section of Jjif/ with transverse momentum 
(/7 t) above 4 GeV/c at mid-rapidity in p -r p collisions at = 500 GeV by the STAR experiment. The ratio of the 
yield of 0^(25) to //0 integrated over A < pj < \2 GeV/c is also presented. Furthermore, the //0 yields are studied in 
different event multiplicity bins in different //0 p-y regions, where the low p-y measurement is enabled by the newly 
installed Muon Telescope Detector. A strong increase of the relative //0 yield with the event multiplicity is observed 
for all /7 t with significant py; dependence. 
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1. Introduction 

In relativistic heavy-ion (HI) collisions, a new state 
of matter, usually referred to as “Quark Gluon Plasma 
(QGP)”, is produced, where quarks and gluons instead 
of hadrons are the relevant degrees of freedom. Quarko- 
nia have been an important probe to study the prop¬ 
erties of QGP as their constituent quarks are predom¬ 
inantly produced in the very early stage of the colli¬ 
sion. In particular, //0, a cc bound state, is predicted 
to be suppressed due to the color screening of its bind¬ 
ing constituent quarks by the surrounding partons in 
the medium, which is viewed as a signature of decon¬ 
finement m. Interpretation of the //0 suppression ob¬ 
served in HI collisions requires thorough understanding 
of its production mechanism in elementary p + p col¬ 
lisions. Despite decades of continuous efforts, under¬ 
standing of the production mechanism still remains in¬ 
complete. Measurements of the //0 production, includ¬ 
ing both the inclusive cross section and its dependence 
on the event multiplicity, at different collision energies 


could shed light on the understanding of different mech¬ 
anisms. In particular, the event multiplicity dependence 
explores the interplay between soft and hard production 
processes, and could provide new insights into the prop¬ 
erties of high-multiplicity p + p events. 

2. Detector and data sample 

The Solenoidal Tracker At RHIC (STAR) 13 is a 
dedicated detector to study HI collisions with coverage 
over full azimuthal angle and \r]\ < 1 in pseudo-rapidity. 
The Time Projection Chamber measures momenta of 
traversing charged particles with good precision down 
to low pj. It also provides ionization energy loss mea¬ 
surement used to identify electrons and muons from //0 
decays. Hadrons can be further distinguished from elec¬ 
trons at high /7 t with the help of the Barrel ElectroMag- 
netic Calorimeter (BEMC), in which electrons deposit 
almost all their energies while the hadrons leave mostly 
only a minimal amount. On the other hand, muon pu¬ 
rity can be further improved using the Muon Telescope 
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Detector (MTD) based on its precise timing measure¬ 
ment as hadrons need more time to reach MTD. Also, 
the magnet in front of the MTD serves as a hadron ab¬ 
sorber. Due to physical constraints, the full MTD covers 
about 45% in azimuth within \ri\ < 0.5. Only 63% of the 
MTD trays were installed in 2013 when the data used in 
this analysis were taken. 

The analysis of reconstructing Jlif/ via the di-electron 
channel utilizes the data sample recorded in 2011 for 
p -h p collisions at = 500 GeV. A High Tower (HT) 
trigger in the BEMC, which requires at least one tower 
with energy deposition above 3.5 GeV, was used to en¬ 
hance the statistics for high It sampled an in¬ 

tegrated luminosity of 22 pb“^ In 2013, the MTD was 
commissioned to trigger on muons from heavy flavor 
decays. Speciflcally for quarkonium measurements, a 
dimuon trigger was used which requires at least two 
MTD hits in coincidence with a MB trigger, and it sam¬ 
pled an integrated luminosity of 28.3 pb“^ data for p -r p 
collisions at 500 GeV. 

3. Results 

3.1. Signal extraction 



Figure 1: Invariant mass distribution of unlike-sign (black histogram) 
and like-sign (shaded area) pairs of electron and position candidates. 
The signal distribution (solid circle) is obtained by subtracting like- 
sign distribution from the unlike-sign. The [1/{2S ) signal is shown in 
the insert as well. 

The invariant mass distribution of unlike-sign 
electron-positron pairs (black histogram) above 4 
GeV/c is shown in Figj^ along with a similar distri¬ 
bution for like-sign electron or positron paris (shaded 
area), which represents a data-driven estimation of the 
combinatorial background. The signal distribution (red 
circles) is obtained by subtracting the like-sign distribu¬ 
tion from the unlike-sign, and is fitted with a Crystal- 
ball function describing the Jlij/ signal shape as well 


as an exponential function describing the residual back¬ 
ground mainly from cc, bb decays and the Drell-Yan 
process. The raw Jlij/ yield is obtained by subtract¬ 
ing the integral of the fitted residual background from 
the sum of the bin contents of the signal distribution 
in the mass range 2.7 < Mee < 3.3 GeV/c^. Loss of 
Jlij/ counts below the accounted mass regime due to the 
missing photon in Jlij/ e^e'y decay is estimated to 
be ~ 10% using the fitted Crystal-ball function. A sim¬ 
ilar procedure is used to extract the raw yield of ij/ilS ) 
in the mass range 3.5 < Mee < 3.8 GeV/c^. 

The invariant mass distribution of unlike-sign muon 
pairs is shown in Fig. (red circles) using the MTD- 
triggered data. Contrary to the BEMC trigger, the MTD 
is able to trigger on low pj muons, which opens access 
to measurements down to zero Jlif/ pj. The invariant 



Figure 2: Invariant mass distribution of unlike-sign muon pairs (red 
circles). It is fitted with a Gaussian function describing signal and an 
exponential function describing background (solid line). 

mass distribution is fitted with a Gaussian function de¬ 
scribing signal and an exponential function describing 
background. The raw Jlij/ yield is extracted by sum¬ 
ming the counts of the unlike-sign distribution in the 
mass range 2.8 < < 3.3 GeV/c^ after removing the 

background contribution estimated via the fitted func¬ 
tion. 

3.2. Inclusive cross section 

Figureshows the inclusive differential cross section 
for Jlij/ multiplied by the branching ratio as a function 
of /7 t at mid-rapidity (|y| < 1) measured through the 
di-electron channel. The total systematic uncertainty, 
shown as open boxes around the data points, is about 
6.3%. It includes two main contributions: i) variation 
of kinematic and electron identification cuts; ii) differ¬ 
ent mass window and different fitting process to extract 
signal. A Non-Relativistic Quantum ChromoDynamics 
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(GeV/c) 

Figure 3: Inclusive Jl\f/ cross section as a function of px at mid¬ 
rapidity (lyl < 1) for p -I- p collisions at 500 GeV measured by STAR 
(red circles). A NRQCD calculation of prompt //j/^ is shown in grey 
band, and agrees reasonably well with data. 


(NRQCD) calculation of prompt JIII/ O is shown as 
the grey band in Fig. Given that only up to 20% of 
the inclusive Jlij/ comes from B-hadron decay HI, the 
agreement between the theoretical calculation and the 
measurement is quite good. 

To help quantify the feed-down contribution of ex¬ 
cited charmonium states, the ratio of the integrated yield 
of 0^(25) within 4 < px < 12 GeV/c to that of //0 
is measured and shown in Fig. It is about 3.5% at 
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Figure 4: Ratio of integrated yield between il/{2S) and Jlif/ within 
4 < pt < 12 GeV/c, \y\ < 1. Similar measurements from other 
experiments are shown for comparison. 


Px ~ 6.8 GeV/c with a total systematic uncertainty of 
about 18%. Compared with other measurements at dif¬ 
ferent colliding energies for p -f p or p -f A systems, 
the current measurement follows the trend of de¬ 
pendence and shows no colliding energy dependence 

ElEllTl. 


3.3. Event activity dependence 

Minimum-bias (MB) p-hp collisions can be classified 
into different event activity bins based on the charged- 
particle multiplicity of the event. To eliminate pile-up 
effects, charged particles are required to leave a signal 
in the Barrel Time-Of-Flight (TOF) detector, which is a 
fast detector insensitive to pile-up. 

To reduce systematic uncertainties and facilitate com¬ 
parison to other experiments and to theory, the TOF- 
matched particle multiplicities and the invariant //0 
yields in different event activity bins are normalized re¬ 
spectively by the same quantity in the MB sample. The 
resulting dependence of the relative //0 yield on the 
event activity is shown in Fig. [^in three //0 pj bins. 
The relative //0 yield is seen to rise monotonically with 
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Figure 5: Dependence of relative invariant JliJj yield on the event ac¬ 
tivity in three different Jjif/ px bins. The boxes around the data points 
represent the systematic uncertainties. 

the increasing event multiplicity, and the rising trend 
is steeper for high px For events with multiplic¬ 
ity about 3 times higher than the average, the relative 
//0 yield above 8 GeV/c is enhanced by a factor of 
about 7.5. For /^x dependence, the relative //0 yield of 
4 < /?x < 8 GeV/c bin is about 62% higher than that of 
/7x > 8 GeV/c bin in the lowest event activity bin, while 
it is about 34% lower in the highest event activity bin. 
The systematic uncertainty is ~ 11% for /^x >0 GeV/c 
bin from the di-muon channel, and ~ 2% for 4 < /7x < 8 
GeV/c and for px > 8 GeV/c bins from the di-electron 
channel. There is an additional one-sided 15% uncer¬ 
tainty which is due to the lack of reliable simulation to 
correct for event selection bias in the MB sample. 

PYTHIA8 |l8l| predictions (v8.183, default tune) for 
//0 pt: above 0 and 4 GeV/c are shown as the dashed 
lines in Fig. with the shaded areas representing the 
statistical errors. They agree quite well with data for 
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both pj bins. The percolation model ii also agrees well 
with data for > 0 GeV/c. However, it is worth not¬ 
ing that the percolation model differs significantly from 
PYTHIA at higher event activity bins, which calls for 
more measurements. 


- p+p collisions @ 500 GeV 
★ STAR: lyl<0.5, p^>0 GeV/c 

^ STAR: |y|<1, p^>4 GeV/c 

^ PYTHIA8.183 default: p^>0 GeV/c 
. PYTHIA8.183 default: p^>4 GeV/c 
— Percolation model: p >0 GeV/c 


STAR preliminary 
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4. Summary 

Differential inclusive Jlif/ cross-section for transverse 
momentum above 4 GeV/c is measured at mid-rapidity 
in p -r p collisions at = 500 GeV via the di-electron 
channel. A NRQCD calculation of prompt Jlij/ produc¬ 
tion agrees well with data in the kinematic range where 
the model is applicable. The ratio of yields of 0^(25) 
to //0 is seen to follow the global trend of /7 t depen¬ 
dence. Furthermore, the relative //0 yield is studied as 
a function of event activity in different py bins. A strong 
correlation between //0 yield and event activity is ob¬ 
served with significant px dependence. Both PYTHIA 
and the percolation model can reproduce the data quite 
well, and the rising trend at RHIC is consistent with that 
observed at the LHC. More measurements at RHIC are 
crucial to distinguish different models, and help under¬ 
stand the characteristics of high multiplicity p-hp events. 


Figure 6: Comparison of relative Jlif/ yield vs. event activity between 
STAR measurement, PYTHIAS prediction and the percolation model. 


Similar measurements for both Jjif/ and D meson are 
observed by the ALICE collaboration for p-f p collisions 
atVTeVlUOlIIIl, shown as open and closed circles in 
Fig. [T] The stronger-than-linear rising trend at RHIC 


- p+p collisions 

- ^ STAR 500 GeV: lyl<0.5, p^>0 GeV/c 

^ STAR 500 GeV: J/v|/^e*e-, |y|<1, p^>4 GeV/c 
O ALICE 7 TeV: |y|<0.9, p^>0 GeV/c 

• ALICE 7 TeV: D meson, |y|<0.5, 2<p^<4 GeV/c 
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Figure 7: Comparison of relative Jlif/ yield vs. event activity between 
STAR (500 GeV) and ALICE measurements (7 TeV). A similar mea¬ 
surement of D meson at 7 TeV by ALICE is also shown. 


follows the trend at the LHC quite nicely, which indi¬ 
cates a possible universal dependence of the relative //0 
yield on the event activity at different energies. More 
differential measurements at different colliding energies 
are needed to verify this observation. 
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